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Abstract 

The relationship between the presence of defects at the stacking structure of the Si/Si02 interface 
and leakage current is theoretically studied by first-principles calculation. I found that the leakage 
current through the interface with dangling bonds is 530 times larger than that without any defects, 
which is expected to lead to dielectric breakdown. The direction of the dangling bonds is closely 
related to the performance of the oxide as an insulator. In addition, it is proved that the termination 
of the dangling bonds by hydrogen atoms is effective for reducing the leakage current. 
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I. INTRODUCTION 



The relevant length scale in silicon-based electronic devices has been continuously de- 
creasing and is currently few nanometers. The understanding of the structural properties 
of Si/Si02 interfaces on the atomic scale acquires an enhanced technological importance 
because the defects at Si/Si02 interfaces cause serious problems, such as the increased_gate 
leakage current, the reduced threshold for dielectric breakdown, and oxide charging M,^- 
These problems are known to degrade the dielectric reliability and device performance. Si02 
films on a Si(OOl) substrate are known to be mostly amorphous. However, intensive and 
sustained scientific activities by various methods up to now have shown that a crystalline 
phase of the Si02 can be observed down to ~ 10 A from the Si/Si02 interface in some 
cases js, [4, 0, Q]. Although the dangling bonds in the amorphous structure of the Si02 
films exhibit the complex and certain effect on the leakage, the influence of the defect at the 
interface becomes more important because atomic-layer Si02 is necessary even under high-/c 
dielectrics. 

The ideal stacking structure of Si/Si02 has an interface consisting of only silicon and 
oxygen atoms without any defects; however this structure is extremely difficult to be formed 
experimentally in the entire interface region. Therefore, the termination of Si-dangling 
bonds associated with O vacancies by some proper atoms is a very important technology. 
Hydrogen has been known to play a crucial role in the fabrication of high-quality Si/Si02 
interfaces. However, the degradation of devices has been attributed to hydrogen releasing 
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from dangling bonds by hot electrons and diffusing to the interface [3, 18| . When the Si02 
layer becomes very thin (~ 15 A), these dangling bonds are expected to increase markedly 
the leakage current through thin oxides as well as to degrade the dielectric reliability due 
to the direct tunneling of electrons from the substrate to the gate electrode. Therefore, 
the relationship between the presence of the defects at the interfaces and leakage current is 
definitely of interest. 

In this paper, a first-principles study on the leakage current of ultrathin gate oxides with 
a Si02 thickness of ~ 14 A is carried out. In many cases, the leakage current thorough 
the Si/Si02 interfaces has been discussed by the one-dimensional Wentzel-Kramers-Brillouin 
approximation and/or the method using empirical parameters, which do not have enough 
accuracy to examine the effect of the dangling bonds from an atomistic point of view. 
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There are also a considerable number of the first-principles studies exploring the dielectric 
breakdown of Si02 films ^] , however the conclusive picture of the relationship between the 
leakage current and the presence of Si-dangling bonds at Si/Si02 interfaces is still missing. 
I found that the resultant leakage current of the interface with the dangling bonds is 530 
times larger than that without the dangling bonds and the length of the charge distribution 
of the defect state in the direction perpendicular to the interface is related to the leakage 
current. In addition, the termination of Si-dangling bonds by hydrogen atoms is effective 
for reducing the leakage current. 

The rest of this paper is organized as follows. In Sec. II, I briefly describe the compu- 
tational methods and models used in this study. My results are presented and discussed in 
Sec. III. I summarize my findings in Sec. IV. 



II. COMPUTATIONAL METHODS AND MODELS 



Our first-principles calculation method for obtaining the electron-conduction properties is 



based on the real-space finite-difference approach 
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12j, which enables us to determine 



the self-consistent electronic ground state with a high degree of accuracy by a timesaving 



double-grid technique [ll|, ll2| and the direct minimization of the energy functional |13| . 



Moreover, the real-space calculations eliminate the serious drawbacks of the conventional 
plane-wave approach, such as its inability to describe nonperiodic systems accurately. I 
examine the electron-conduction properties of the Si/Si02 interface suspended between semi- 
infinite electrodes by computing the scattering wave functions continuing from one electrode 
to the other. 

Figure [1] shows an example of the computational model for electron-conduction calcula- 
tion. The Si/Si02 interface structure proposed by Kageshima and Shiraishi |l^ in which 
an a-quartz crystalline Si02 layer stacked on a Si(OOl) surface is employed. From the first- 
principles analysis for the interface energy, the interface structure employed here is one of 
the candidates among the possible interface structures ISj. The model shown in Fig. [1] cor- 
responds to an interface without any dangling bonds, which I refer to as the perfect interface 
hereafter. The thickness of the Si02 layers is ~ 14 A and that of the silicon substrate is 
7.7 A. According to the oxidation process proposed by Kageshima and Shiraishi IJ], the 
perfect interface appears after the emission of a silicon atom when six oxygen atoms are in- 
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serted into a silicon bulk. However, even when the number of inserted oxygen atoms is five, 
the silicon atom can be emitted and two dangling bonds remain. Since hydrogen is known 
to passivate the dangling bonds, I employ the interface models assuming that one oxygen 
atom at the interface is absent and the dangling bonds are terminated by hydrogen atoms. 
Figure [2] shows the computational models: the perfect interface without any defects [model 
(a)], the hydrogen-atom bridge parallel to the interface [model (b)], Pb center [model (c)], 
the hydrogen-atom bridge perpendicular to the interface [model (d)], the hydrogen- molecule 
bridge parallel to the interface [model (e)], and the hydrogen-molecule bridge perpendicular 
to the interface [model (f)]. Model (b) [models (c) and (d)] is expected to emerge after one 
hydrogen atom in model (e) [model (f)] is released from dangling bonds by hot electrons. 
The side length of the supercell is taken to be 10.86 Ax 5.43 A and A;-space integrations 
are performed with 1 x 2k points in the irreducible wedge of a two-dimensional Brillouin 
zone. The norm-conserving pseudopotentials 17| of Troullier and Martins are employed 
to describe the electron-ion interaction. Exchange correlation effects are treated by local 
density approximation lOj and freedom of spin is not considered. 



III. RESULTS AND DISCUSSION 



To determine the atomic configuration of the interface, I first perform structural optimiza- 
tion using a thin film model in which the dangling bonds of silicon atoms in the bottommost 
and topmost layers are terminated by hydrogen atoms [2^. The numbers of atoms in the 
models are listed in Table [H The periodic boundary condition is imposed in the x and y di- 
rections, while the isolated boundary condition is employed in the z direction. The vacuum 
region above and beneath the thin film is ~ 6.5 A. I relax all the atoms in model (a) to 
determine the thickness of the thin film and then execute structural optimizations for the 
other models. For models (b), (c), (d), (e) and (f), the Si atoms in the first bottommost, 
second bottommost, and topmost layers, and the H atoms passivating the dangling bonds 
of the Si atoms in both the surface layers are frozen at the positions obtained by model 
(a) during structural optimization. Using the above parameters, formation energies of the 
interface defects are calculated for their optimized geometries. The fifth column of Table [T] 
shows the formation energies Ef defined as 

Ef = + (33 - No)fio + {Nh - l2)fiH - E^^)- (1) 
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Here, -E'(a), -E'(x), A^o, and hh are the total energies of model (a) and the other models, and 
the chemical potentials of oxygen and hydrogen atoms, respectively. In addition. No (Nh) 
is the number of oxygen (hydrogen) atoms. The chemical potentials of oxygen and hydrogen 
atoms are set to be those of their dimers, however the definition of them does not affect 
my conclusion as long as I compare the energies of the models with the same numbers of 
No and Nh- The atomic configuration of model (d) is energetically metastable and the 
total energy of model (d) is higher than those of models (b) and (c) by ~ 0.6 eV/supercell. 
In addition, one can find that the hydrogen atoms can be released easily as expected from 
experiments by comparing the formation energy between the models with and without the 
dangling bond, although the quantitative evaluation is difficult because the bonding energy 
is affected by the definition of the chemical potential. 

Next, to examine the electron-conduction properties, I remove the H atoms terminating 
the dangling bonds of both the surface and suspend the Si/Si02 interface between aluminum 
jellium electrodes, as shown in Fig. [H The distance between the Si atoms and the edge of 
the jellium electrodes is chosen to be 2.42 A so as to imitate the Si(OOl) and Al(lll) 
interface. Figure [3] shows the charge-density distribution of the states between Ep—O.b eV 
and Ep+0.5 eV. There are no gap states in models (a), (e), and (f), while the states due to 
dangling bonds are observed in models (b), (c), and (d). The charge-density distributions of 



the gap states are spatially localized around the defects and agree wit 
occupied orbital of the neutral hydrogen bridge [models (b) and (d)] 



]__that of the highest 
and with that of 



the highest occupied orbital of the Pb center [model (c)] [22] . In model (c), the surface silicon 
atom having a dangling bond moves down to the silicon substrate so that its bond angles 
become planar and the vr electron forms a Pb center in the silicon substrate. On the other 
hand, the dangling bond in model (d) forms vertically and spreads into the dioxide. The 
electron-conduction properties of the interface system are computed using the scattering 
wave functions continuing from one electrode to the other, which are evaluated by the 



overbridging boundary- matching method Uj, l23|| under the semi- infinite boundary condition 



in the z direction. The leakage current of the system at the limits of zero temperature is 



evaluated using the Landauer-Biittiker formula 



/ = --y^ nE)dE, (2) 

where T{E) is the transmission probability of the incident electron from the lower electrode 



with the energy E. The leakage current is computed by the integration of the tunnel current 
over the range from Ei = Ep — 0.5 eV to E2 = Ep + 0.5 eV, where Ep is the Fermi energy 



25|. 



The sixth column of Table [I] shows the leakage current ratio in which the current of the 
perfect interface is set to be unity [26|. Note that the passivation of the dangling bonds 
arising from the oxygen vacancy by the hydrogen molecule markedly reduces the leakage 
current. On the other hand, the current increases with the presence of dangling bonds at 
the interface. In particular, the current in model (d) is 534.3 times larger than that in model 
(a) and the defect in model (d) is expected to cause serious problems in the practical devices. 

To explore the relationship between the presence of dangling bonds and leakage current 
in detail, I plot in Fig. H] the charge- density distribution of the scattering wave functions 
with the Fermi energy. Most incident electrons from the Si substrate are reflected at the 
Si/Si02 interface. In addition, the characteristics of the charge-density distribution around 
the defects completely agree with those of the gap state in Figs. [3|b), [3](c), and[3](d). Thus, 
the incident electrons penetrate through the gap states. On the other hand, I could not find 
the evidence that the dangling bonds terminated by the hydrogen atoms cause the leakage 
although the charge- density distributions of the scattering wave functions look similar with 
those in Figs. [3](e) and[3](f). The dangling bond states energetically shift up to the valence 
band or down to the conduction band by the passivation so that the states do not contribute 
to the leakage. Figure [5] shows the densities of states (DOS) of models (a), (b), (c) and 
(d), which are plotted by integrating them on the plane parallel to the interface following 
p{z,E) = J \'ilj{r,E)\'^dr\\. The width of the energy band gap at the Fermi level around 
the interface region in model (a) agrees with that obtained by the other first-principles 



calculation, in which more thick Si and Si02 films are directly bonded [27|]. This implies 
that thickness of both the film in my model is sufficient to eliminate unfavorable effect 
from the electrode. Most electrons of the gap states in model (b) accumulate beneath the 
interface atomic layer because the direction of the dangling bond is parallel to the interface. 
In Fig. [5](c), the two peaks in the distribution of the DOS at the Fermi level is the evidence 
that the gap state of model (c) comes from the Pf, center, which is composed of the tt 
electrons. However, the state locates spatially deep inside the silicon substrate and does not 
contribute to the leakage significantly. In contrast, the charge-density distribution of the gap 
state in model (d) spreads along the direction perpendicular to the interface more widely 
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than those of the other models and penetrates into the oxide region because the hydrogen- 
atom bridge ahgns perpendicular to the interface and the silicon atom having dangling bond 
exists inside the oxide. According to the problem of the tunnel effect of a potential barrier in 
quantum mechanics, the transmission probability exponentially decreases with an increase 
in the potential-barrier length. The dangling bond in model (d), which locates in the energy 
band gap of the Si02 film and spreads in the direction perpendicular to the interface, results 
in a significant increase in leakage current. Thus, it degrades the performance of the oxide 
film as an insulator film of metal oxide semiconductor field-effect transistors. 

IV. CONCLUSION 

1 have implemented first-principles simulations to examine the relationship between the 
defects at the Si/Si02 interface and leakage current. I found that the Pb center and the 
hydrogen- atom bridge aligned parallel to the interface do not affect the leakage current 
significantly. On the other hand, because of the spread of the charge-density distribution 
of the gap state, the hydrogen-atom bridge aligned perpendicular to the interface results 
in the largest leakage current, and therefore, it causes dielectric breakdown. In addition, 
the hydrogen termination of the dangling bonds significantly reduces decrease the leakage 
current. 
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TABLE L Number of atoms in supercell, formation energies and leakage current ratio. The models 
correspond to those shown in Fig. [2j 



Model 


Si H 


Ef (eV) 


Current ratio 


(a) 


38 33 12 




1 


(b) 


38 32 13 


3.66 


3.9 


(c) 


38 32 13 


3.59 


12.6 


(d) 


38 32 13 


4.23 


534.3 


(e) 


38 32 14 


3.21 


0.9 


(f) 


38 32 14 


3.01 


1.6 
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FIG. 1: Schematic image of computational model. Gray and black circles represent silicon and 
oxygen atoms, respectively. The rectangle enclosed by broken lines represents the scattering region 
when the leakage current is calculated. 




FIG. 2: Si/Si02 interface models. Perfect interface without any defects (a), hydrogen-atom 
bridge parallel to interface (b), Pb center (c), hydrogen-atom bridge perpendicular to interface 
(d), hydrogen- molecule bridge parallel to interface (e), and hydrogen-molecule bridge perpendicu- 
lar to interface (f). The meanings of the symbols are the same as those in Fig. [T]and white circles 
represent hydrogen atoms. 
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FIG. 3: Contour plots of charge-density distributions for electrons with energies between Ep — 0.5 
eV and Ef+0.5 eV. The planes shown are along the cross section in the (110) plane of the silicon 
substrate including dangling bonds. Each contour represents twice or half the density of adjacent 
contour lines and the lowest contour is 4.22 x 10"*^ e/A^. The meanings of the symbols are the 
same as those in Fig. [2j (a), (b), (c), (d), (e), and (f) correspond to models (a), (b), (c), (d), (e), 
and (f ) in Fig. [2l respectively. 
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FIG. 4: (color) Charge-density distributions of scattering wave functions with the Fermi energy. 
Each contour represents twice or half the density of adjacent contour lines and the lowest contour 
is 4.72 xlQ-^ e/eV/A^. (a), (b), (c), (d), (e), and (f) correspond to models (a), (b), (c), (d), 
(e), and (f) in Fig. [21 respectively. Light blue, red, and white balls represent silicon, oxygen, and 
hydrogen atoms, respectively. 
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FIG. 5: Distributions of DOS integrated on plane parallel to interface as functions of relative energy 
from the Fermi energy. Zero energy is chosen to be the Fermi energy. Each contour represents 
twice or half the density of adjacent contour lines and the lowest contour is 1.39 xlO~^ e/eV/A. 
(a), (b), (c), and (d) correspond to models (a), (b), (c), and (d) in Fig. [21 respectively. Dashed 
and dotted lines represent the vertical positions of silicon and oxygen atomic layers, respectively. 
Arrows denote the oxygen atomic layer from which an oxygen atom is extracted. 
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